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1 Introduction 



B 

■ Quantum Hall ferromagnets at filling factor u = 1 exhibit charged excitations that 

Ch , are topologically stable objects called skyrmions. Various experiments [jl], |^ have re- 

Q [ vealed strong evidence for the existence of these objects in two-dimensional electron 

systems. 

For their theoretical description we employ the microscopic Hartree-Fock (HF) 
theory with the hedgehog ansatz for charged spin-textures, developed by Fertig et. 
al. [P] and subsequently commonly used in the literature |^, We show that the 
most general single Slater determinant in the lowest Landau level (LLL) approxima- 
tion collapses into the hedgehog ansatz, hence supporting this approach. Physical 
observables like magnetization and charge density are then calculated microscopi- 
cally in the standard way (section 2). However in contrast to previous approaches, 
we will analyse the full magnetization vector, not only its third component, the spin- 
density. It will turn out that the magnetization is not a unit vector, but its length 
is substantially decreased in regions with non-zero excess charge (bag formation). 
A signature of this behaviour was already observed in 0, namely that the absolute 
value of the spin-density at the origin as obtained from the HF calculation is reduced 
compared to that of the non-linear 0(3) model which is always unity. However the 
conjecture that with quantum corrections included, the discrepancies between ef- 
fective field theory and microscopic HF results may be resolved, is incorrect. Bag 
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formation appears already in tree level and accordingly the magnetization vector 
can no longer be identified directly with the 0(3) unit vector field. An additional 
scalar field is necessary to describe the magnetization vector's length. 

Subsequently, using the gradient expansion, a new 0(3) type of effective action 
for extended spin-textures is derived which takes these effects into account (section 
3) . It will be demonstrated, that this action reproduces the results of the microscopic 
HF calculation correctly for zero range forces. 

2 Microscopic Hartree-Fock calculations 

As usual lengths are measured in units of the magnetic length £ = ^UcjeB and 
energies in units of the Coulomb energy e^/el. Then for a given electron-electron 
interaction in the LLL there exists only one dimensionless parameter in the theory, 
namely the effective Zeeman coupling g = Qs^bB / {e^ / et). 

In the completely filled LLL (filling factor v = 1) the ground state is given by 
the unique single Slater determinant (Vandermonde's determinant, see e.g. 0) 

oo 11 

1^) = ^ n ^ ^r\m) = ^^,r-e--e-^'/' (1) 

m=0 ^ 

with the spins of all single particle states down. Here we are concerned with the de- 
scription of charged spin-texture excitations, in particular with those in the charge 
sectors Q = ±1. Due to the particle-hole symmetry the energies of the posi- 
tively charged (quasi-hole, antiskyrmion) and negatively charged (quasi-particle and 
skyrmion) excitations, = E^eo, differ only by the exchange energy per electron 

in the ground state, e.g. eo = —\Jtx/S in the pure Coulomb case. Thus it is sufficient 
to consider the Q = +1 sector with the following ansatze for the many-body wave 
function: 

(i) the quasi-hole 

oo 

1^) = n l'^ + 1 i) ' (quasi-hole) (2) 

m=0 

with just one electron removed from the completely filled LLL. Its energy relative 
to the ground state is E^ = g/2 — 2eo {E = g/2 — eo). 

(ii) the hedgehog 

oo 

1^) =^ n (^-l"^t)+^^m|"^ + l i)) , ul + vl = l, (hedgehog) (3) 

m=0 

introduced by Fertig et. al. 0. For Um = 0, Vm = Vm, this ansatz reduces to the 
quasi-hole (0). 
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(iii) the general single Slater determinant 

\^)= Allies), |«) = E C:,Jma) , (general single SD) (4) 

a ma 

whose single particle states \a) with {a\a) = 1 are given by the most general linear 
combination of LLL-states labeled by magnetic quantum number m = 0, . . . , N — 1 
and spin a (up and down). 

A HF calculation for the general Slater determinant (ansatz (|^)) was performed 
with — 1 particles in order to guarantee Q = +1 in a basis up to = 300 states 
involving the diagonalization of 600 x 600 matrices. Independently from the initial 
conditions chosen, the general Slater determinant always collapsed 

^ml ~ '^a^m,a ; ^mj ~ 'Va^m,a+1 (5) 

into the specific hedgehog ansatz (^ during the iteration process. This happened for 
various Zeeman coupling strengths g and also for different (repulsive) interactions. 
Thus, we may consider this finding as a strong indication that the lowest energy 
single Slater determinant for Q = 1 in the LLL is indeed always of the hedgehog 
form. I am not aware of a rigorous proof of that statement, although it is certainly 
related to the fact that in the ansatz (|^) only single particle states with the same 
grand-spin m + a couple. 

For the reasons just explained we may use the hedgehog ansatz (|^) in the Q = 
+ 1 sector which simplifies the calculation considerably. Instead of diagonalizing 
matrices there are only coupled algebraic equations for Um and Vm to solve for each 
m separately. We give the expression for the HF energy derived in |Q relative to 
the completely filled LLL 

E^F = fE(^^-^m-i + l) 

m 

+ ^ E("m + ^m-l - l)(«m' +^m'-l " I) {m TTl' \V\m m) 
mm' 

I 

mm' 

— UmVmUm'Vm' fn' + l\V\m' m + 1) . 

mm' 

The individual terms are the Zeeman energy and the direct and exchange terms 
of the two-body interaction V. The corresponding HF equations are immediately 
derived. The (excess) charge density 

p^(r) ^ ^-{^\j:^{r-rm 

1 2 
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is radially symmetric and normalized to Q = +1 as may be checked using the 
condition + t>^ = 1. The magnetization vector 

m(r) = — 27r(^| ^ cr(i)5(r — Ti)|^E') = [^fij) cosy9, /(r) sin;/), /i(r)) 

/(r) = -V2rE^^^e-% = E(^™-i " ^ (§) 

^ y/m + 1 m! ^ m! 

{cr{i) are the Pauli matrices for particle i) normalized to lim^.^oo'^T' = ^3 is of the 
hedgehog form as expected. However, it is certainly not a unit vector: there are 
two radial functions involved for which in general + /i^ 7^ 1. In order to make 
closer contact to effective field theories we may rewrite the magnetization vector in 
the form 

m{r) = o"(r)Th = cr{r) (sin F{r) cos (p, sin F{r) sin cp, cos F{r)) (9) 

with rh a unit vector. The scalar function a describes the deviation of the magneti- 
zation vector's length from 1 (bag formation) and the angle function F determines 
the orientation of the unit vector rh. It is clear that these functions are readily 
expressed in terms of the original functions / and h. 

In this way we are going to analyse the full magnetization vector m, not only 
its third component, which is the main difference to previous presentations. It will 
immediately lead us to an important relation of the formed bag with the charge. In 
order to better understand the situation we will study the analytical hard-core model 
before we present numerical results for the more realistic Coulomb interaction. 



2.1 Hard-core interaction 

For a zero range interaction V(ri — 7^2) = t5(ri — 7^2) and vanishing Zeeman splitting, 
g = 0, the HF equations may be solved analytically 



I a / m + 1 , , 

^rn = \ , , , , Vrr, = \ , a=—. 10 

ym+l + a ^ m + 1 + a 2 

Here, A is a free parameter which fixes the soliton size. The radial functions deter- 
mining the magnetization m in (^) are 

f{r) = —\rp{x,a), h{r) = {x — a) p{x,a) , x = — (11) 
with the single valued analytical function 

1 x^ 
p{x,a) = J2 — r^. re"'', 

x"+^p(x,a)e^ = rdtf'e', (12) 



— {xp{x, a)) + {a + x)p{x, a) 
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related to the incomplete gamma function . From and (0) we obtain 



F(r) = 2 arctan — , a(r) = {x + a) p{x, a) (13) 
r 

the angle function of the familiar Belavin-Polyakov soliton [Q, but the modulus 
function a is not equal to 1. The topological and physical charge densities are 

P^(r) = -^-^(cosF) = „ ,„ , (14) 
p''{r) = ^\\-{a\x)v{x,a)\=^—^. (15) 



It is readily checked with (|T^) that both charges are = = +1, however the 
densities do look quite different in general. For example for the quasi-hole, A — > 0, 



we obtain from ([12|) a = \ — exp(— r^/2), i.e. 

p^(r) = 5(r), A^O, (16) 

p''(r-) = ^exp(-rV2), A^O. (17) 



It is only in the opposite limit of very large skyrmions A — > cxd that the two densities 
become similar. But note, that although o \ WxqXj limit, the integral over the 
physical charge density, i.e. over (1 — a)/2Ti, always remains 1. 

Thus, the conclusion of this simple exercise is twofold: (i) p + g'^ = a'^ ^ 1, 
a bag is formed which does go beyond the standard 0(3) description and (ii) the 
physical charge density is not equal to the 0(3) model's topological density, rather 
it is connected with the bag via (plSj) . 

It should be added, that for hard-core forces at finite Zeeman splitting the size 



parameter A and the (anti-) skyrmion collapses into the quasi-hole (|T6|,|T^. 
There are no skyrmions for hard-core forces. In the following we are going to discuss 
the situation for Coulomb forces. 



2.2 Coulomb interaction 

Although a realistic calculation would have to take into account Landau level mixing, 
the finite thickness of the layer 0, [1^, ^, |lT| and also impurities, we consider the 
electrons in the LLL interacting in the plane via pure Coulomb forces V{ri — r2) = 
1/ 1 ri — r2 1 . We are interested in bag formation and its relation to the physical charge 
density. The corresponding HF equations are solved numerically using a basis up 
to several thousand states for large skyrmions. The soliton energy in dependence 
of the Zeeman coupling strength g is depicted in Fig. 1. We observe a second 
order phase-transition at ~ 0.053 below which (anti-) skyrmions are energetically 
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Figure 1: Skyrmion and quasi-hole energies as a function of the Zccman splitting 
calculated in Hartrcc Fock for a pure Coulomb interaction. For g < 0.053 skyrmions 
are energetically favored. In the limit g ^ the Bogomol'nyi bound is approached. 
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Figure 2: Profiles for the (anti-) skyrmion with gi = 0.03 versus the radius in magnetic 
lengths. A bag is formed. 
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Figure 3: Topological and physical charge densities for the soliton with g = 0.03 
(solid lines). Although both total charges are normalized to 1, the physical charge 
density deviates significantly from the topological one. Instead it is extremely well 
approximated by the bag's shape 1 — a (dashed line). 



favored compared to the quasi-hole (dashed line). For extremely small Zeeman 
couplings (large skyrmions) the soliton energy approaches the Bogomol'nyi bound 
~ .3133, indicating that in this limit the Belavin-Polyakov soliton of the 
non-linear sigma model is obtained. 

The profiles are calculated microscopically via the magnetization (^) and plot- 
ted in Fig. 2 for Zeeman coupling g = 0.03. The formation of a bag is clearly 
indicated. With decreasing Zeeman coupling the bag begins to fiatten out. But 
notice, that in the Belavin-Polyakov limit g ^ the integral over (1 — cr)/27r is 1. 
The physical charge density, calculated microscopically according to (|^), and the 
topological density defined via the angle function F (Fig. 2) are depicted in Fig. 3 
again for Zeeman coupling g = 0.03. It is noticed that the physical charge density 
is much fiatter compared to the topological density. The bag shape function 1 — cr 
(dashed line) is hardly distinguishable from the physical charge density. This sup- 
ports our findings for hard-core forces: it is the bag and not the topological charge 
density which is closely related to the physical charge density. 



3 Effective field theory 

There has been a lot of effort to derive an effective field theory for the quantum Hall 
ferromagnet |T2|, |T^, [T^, |r^, . The result in terms of a gradient expansion of the 
unit magnetization vector rfi is essentially an 0(3) model, i.e. the non-linear sigma 
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model augmented by higher order terms stemming from the interaction. However 
from the above considerations it is expected that the full magnetization vector m, 
namely its modulus a, should show up in the effective action. In particular the rela- 
tion between the physical charge density and the bag ([T5|) should be respected. The 
derivation of such an effective action in full generality is a challenging task. There- 
fore, in the following subsection we resort to the gradient expansion for extended 
spin-textures quite as in previous approaches, however we take bag formation into 
account. 

3.1 0{3)- model 

Starting point is the full magnetization vector and the physical charge density which 
in lowest order gradient expansion read 

m = arh = ^1 — '^diifidi'rfij rfi + (9(4) , (18) 

= -^eijrh{d,rh x d^rh) + 0(4) . (19) 

Because of the charge conjugation properties of the magnetization vector, i.e. the 
scalar field is the same for skyrmions and anti-skyrmions, it is clear that in (^) the 
energy density of the non-linear sigma model rather than the topological density 
enters. It is remarkable that its coefficient, which follows e.g. from ([15D , is indepen- 
dent from the specific two-body interaction. From the discussion in section 2 (cf. 
Fig. 3) it is evident that the physical charge density in general is not well approx- 
imated by the topological density, nevertheless in terms of the gradient expansion 
for extended spin-textures eq. ( p^ holds. In accordance with eq.(|l8]) the energy 
functional could be chosen 

E[m] ~ ^ J d'^r[dimdirn + (1 — mm)^] + mteiaction terms 

= ^ j dh[a^dMrn+{l-a^f] + ... , (20) 



with ps representing the spin-stiffness, e.g. Anps = \J'k/?>2 for the Coulomb inter- 
action. Of course other potentials for the scalar field are possible similar to the 
one derived from the QCD scaling behaviour to describe bag formation in hadron 



physics IIT3] • However up to fourth order in the gradient expansion the exact form of 
the potential will not matter provided that it reproduces ([T8|) correctly. Note in this 
context that gradients on the scalar field are of higher order and may be neglected. 

None of the previously proposed effective actions allows for bag formation. The 
collective approaches [|l^, 0, |l^, 0] operate with unitary rotations which restrict 
the magnetization to the unit vector rh right from the beginning. Derivations from 
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Chern-Simons theories have a scalar field included which however is not taken seri- 



ously dynamical field [0. 

In order to arrive at a non-linear 0(3) type of model, the scalar field may now 
be eliminated from (|20|) in favour of a symmetric fourth order term with definite 
strength related to the spin-stiffness 



Ps f 1 

E[rh] = ~^ J d'^r[dirhdirn — -{dirhdirhY + 0{6)] + interaction terms . 



(21) 



It should be emphasized that although the above 0(3) model deals with a unit 



vector m the bag is still formed via the relation to the magnetization 
(|18D and (^1]) are at variance with the corresponding equations in [p!5[ 
energy functional, including interaction and Zeeman terms is given by 



Eqs. 
The full 



E[rh] 



diThdjiri 



+ 



+ 



rfVirfV2p^(ri)\/(ri - r2)p^(r2) 



(22) 



47r 



1 — ( 1 — -dirhdiiri 



with p*^ replaced by (|T9|) in lowest order gradient expansion. This energy func- 
tional is clearly charge conjugation invariant as it should be. It is unique up to 0{4) 
in the gradients and linear terms in the Zeeman coupling. The essential difference to 
previously proposed effective field theories is the additional symmetric fourth order 
term and of course the nontrivial connection of the 0(3) field with the magnetization 



The energy functional (0) captures the essential features of the microscopic 
approach discussed in section 2. For zero range forces, Airps = t/Air, and vanishing 
Zeeman coupling it simplifies with the hedgehog-ansatz (y) inserted 



E\F] = f^ld' 




sin^F 



(23) 



It is immediately verified that the Belavin-Polyakov soliton (|T3D solves the corre- 
sponding stability condition. With this solution inserted the second term in (|23|) 
vanishes leaving only the pure non-linear sigma model. Hence, the soliton energy 
is E = 4:7Tps = i:/47r as in the microscopic calculation. Also, the charge is related to 
the integral over the bag 



2tt 



in 



sin^F' 



(24) 



as in subsection 2.1. Again, for finite Zeeman coupling the soliton collapses. None 
of the previously proposed effective actions has this simple property that it gives 
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the correct result for extended spin-textures with hard-core forces. It is clear that 
in order to obtain these results, the symmetric fourth order term in (^) is essential. 

In the Coulomb case, we obtain also solitons for non-zero but weak Zeeman cou- 
plings. Because the repulsive interaction term scales with the inverse soliton size 1/A 
it competes with the attractive symmetric fourth order term proportional to — 
and there exists a critical soliton size, determined by the Zeeman coupling strength, 
below which the soliton ceases to exist. In this way the effective theory models 
qualitatively the phase-transition discussed in subsection 2.2. However, because of 
the gradient expansion employed, it is not expected that the energy functional (|22|) 
remains valid for too small solitons. Down to which soliton sizes (^2]) may give a 
fair description of the microscopic HF calculation has still to be investigated. 

It the 1-soliton sector it appears that the effective field theory has not much 
advantage over the microscopic approach. The numerical difficulty in solving the 
non-linear (integro-) differential equation for the profile function deduced from (|22|) 
is comparable to that of the HF calculation. However for multi-soliton systems 
such as skyrmion gases and crystals [|TB|, [T^, the virtues of the effective field 



theory become apparent. In terms of collective coordinates the dynamics of such 
multi-soliton systems may efficiently be described. Also beyond tree approximation 
when loop corrections are considered the effective field theory approach might have 
advantages. On the microscopic level this corresponds to go beyond HF which is 
complicated. 

In the following we give a brief argument why properties like bag formation 
discussed in the context of the HF calculation are not explained by loop corrections. 



3.2 Loop corrections 

Intuitively it is clear that the HF calculation on the microscopic level should corre- 
spond to the tree approximation of the effective theory. Accordingly, bag formation 
( p!8D and the symmetric fourth order term in ( pT]) appear already classically. Could 
quantum corrections be responsible for these effects as suggested in 0? 

In the absence of an appropriate expansion parameter like the number of colors 
in hadron physics the role of loop corrections is difficult to judge. In fact, loop 
corrections may lead to bag formation by softening the constraint = 1. Also 
the energy is lowered and the charge distribution broadened. Using heat kernel 
expansion methods one can even show that loop corrections effectively lead to a 
symmetric fourth order term with the correct (attractive) sign. Nevertheless, we 
argue that these effects observed in the HF calculation cannot be attributed to loop 
corrections for the subsequent argument. 

Let us for simplicity consider only the non-linear sigma model part in the effective 
action. The equations for the fluctuations are independent of the strength of that 
term, i.e. the spin-stiffness ps- The only dimensionful parameter which appears is 



10 



the size A of the BP sohton. For that reason the effective symmetric fourth order 
term produced by the heat kernel expansion comes with a coefficient proportional to 
A such that the total term scales like 1/A. In fact, the whole Casimir energy scales 



like 1/A as was explicitely demonstrated in |21|. In any case, the spin-stiffness does 
not appear at all, which implies that the symmetric fourth order term found in (|2T| ) 
necessary to model the microscopic HF calculation cannot be explained by quantum 
fluctuations. From that argument it becomes clear that the bag is already formed 



in tree level according to (|T^) and (^Tj). 

The above reasoning does not imply that fluctuations are unimportant. In fact, 
we expect significant loop corrections which among other things also will in- 



tensify the bag formation. But in the microscopic picture this would correspond to 
go beyond HF |T0[ performing e.g. a RPA calculation. 



4 Conclusions 

Charged skyrmions in quantum Hall ferromagnets at filling factor u = 1 were re- 
analysed in the framework of the HF theory. It was demonstrated, that within the 
LLL approximation the hedgehog ansatz leads to the lowest energy configuration in 
the unit charge sector. Calculating all three components of the magnetization vector 
it was noticed that its length deviates locally from unity, i.e. a bag is formed which 
accommodates the excess charge. The physical charge density is extremely well ap- 
proximated by the bag, but in general deviates considerably from the topological 
density. Hence, the commonly accepted and celebrated equivalence of physical and 
topological charge does not hold locally. 

In terms of a gradient expansion an improved 0(3) type of effective action is 
proposed for the quantum Hall ferromagnet, which takes bag formation into account. 
Accordingly, the magnetization vector is nontrivially related to the 0(3) unit vector 
field. In addition to previous approaches, the effective action contains an attractive 
symmetric fourth order term. For hard-core forces it was explicitely demonstrated 
that with these amendments made the correct microscopic result for extended spin- 
textures is obtained. For more realistic forces the effective action models the phase- 
transition anti-skyrmion to quasi-hole (skyrmion to quasi-particle) observed in the 
microscopic calculation when the Zeeman splitting is increased. The additional 
fourth order term also modifies the soliton-soliton interaction, e.g. it influences 
the dipole-strength in the dipole-dipole force for skyrmions. By this means the 
properties of multi-skyrmion systems such as skyrmion gas and lattice are affected. 

The present investigation can be completed in several respects. The range of va- 
lidity of the proposed energy functional for realistic two-body forces may carefully be 
tested by comparison with the microscopic HF calculation. Also the time-derivative 



part of the effective action |15] is of interest and will be subject to changes due to the 



bag. Beyond that, the derivation of a corresponding effective field theory without 
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resorting to the gradient expansion is desirable. Such a theory should, in addition 
to the 0(3) unit vector field, contain a scalar field for the magnetization vector's 
length as a dynamical field. Finally, loop corrections may be calculated along the 
lines of ref. I^Tj. On the microscopic level this can be complemented by RPA type 
calculations which take fiuctuations into account. 

5 Acknowledgments 

The author wishes to thank his colleagues G. Holzwarth, T. FlieBbach, F. Meier and 
A.-S. Marculescu for contributing to a seminar about the quantum Hall effect from 
which this study has emerged. This work was supported in part by funds provided 
by the Deutsche Forschungsgemeinschaft (Contract No. DFG-Ho-527/16-1). 



References 



[1] 

[2] 

[3] 

[4] 
[5] 

[6] 

[7] 
[8] 

[9] 
[10] 



A. Schmeller, J. P. Eisenstein, L.N. Pfeiffer and K.W. West, Phys. Rev. Lett. 
75, 4290 (1995). 

S.P. Shukla, M. Shayegan, S.R. Parihar, S.A. Lyon, N.R. Cooper, and A. A 
Kiselev, Phys. Rev. B61, 4469 (2000). 

H.A. Fertig, L. Brey, R. Cote, and A.H. MacDonald, Phys. Rev. B50, 11018 
(1994). 

N.R. Cooper, Phys. Rev. B55, 1934 (1997). 

H.A. Fertig, L. Brey, R. Cote, A.H. MacDonald, A. Karlhede, and S.L. Sondhi, 
Phys. Rev. B55, 10671 (1997). 

M. Abolfath, J.J. Palacios, H.A. Fertig, S.M. Girvin and A.H. MacDonald, 
Phys. Rev. B56, 6795 (1997). 



S.M. Girvin, Les Houches lectures, |cond-mat/9907002 



Ed. M. Abramowitz and A. Stegun, Handbook of Mathematical Functions, 
Dover Pub. Inc., New York, 1965. 

A. A. Belavin and A.M. Polyakov, JETP Lett. 22, 245 (1975). 

B. Kralik, A.M. Rappe, and S.G. Louie, Phys. Rev. B52, 11626 (1995), B56, 
4760 (1997). 



[11] V. Melik-Alaverdian, N.E. Bonesteel, and G.Ortiz, Phys. Rev. B60, R8501 
(1999). 



12 



[12] E. Fradkin, Field Theory of Condensed Matter Systems, Addison- Wesley, Red- 
wood City, 1991. 

[13] K. Moon, H. Mori, K. Yang, S.M. Gervin, and A.H. MacDonald, Phys. Rev. 
B51, 5138 (1995). 

[14] Y.A. Bychkov, T. Maniv, and I.D. Vagner, Phys. Rev. B53, 10148 (1996). 
[15] W. Apel and Y.A. Bychkov, Phys. Rev. Lett. 78, 2188 (1997). 
[16] K. Moon and K. Mullen, Phys. Rev. B57, 14833 (1998). 

[17] H. Gomm, P. Jain, R. Johnson, and J. Schechter, Phys. Rev. D33, 3476 (1986). 

[18] R. Cote, A.H. MacDonald, L. Brey, H.A. Fertig, S.M. Girvin, and H.T.C. Stoof, 
Phys. Rev. Lett. 78, 4825 (1997). 

[19] Yu.V. Nazarov and A.V. Khaetskii,, Phys. Rev. Lett. 80, 576 (1998). 

[20] K. Moon and K. Mullen, |cond-mat/9908310| . 



[21] H. Walliser and G. Holzwarth, Phys. Rev. B61, 2819 (2000). 



13 



